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We extend the saturation models a la Golec-Biernat and Wiisthoff to cross-sections of hard pro- 
cesses initiated by virtual-gluon probes separated by large rapidity intervals at hadron colliders. We 
derive their analytic expressions and apply them to physical examples, such as saturation effects 
for Mueller-Navelet forward jets. By comparison to 7*— 7* cross-sections we find a more abrupt 
transition to saturation. We propose to study observables with a potentially clear saturation signal 
and to use heavy vector and fiavored mesons as alternative hard probes to forward jets. 
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The saturation regime describes the high-density phase of partons in perturbative QCD. It may occur for instance 
- - - when the Balitsky-Fadin-Kuraev-Lipatov (BFKL) QCD evolution equation Q goes beyond some energy limit 2, 

\ 0, 0) 0) • On a phenomenological ground, a well-known saturation model Q by Golec-Biernat and WiisthofF (GBW) 
3^ ■ gives a parametrisation of the proton structure functions already in the HERA energy range. It provides a simple 
' and elegant formulation of the transition to saturation. However, there does not yet exist a clear confirmation of 
^ • saturation since the same data can well be explained within the conventional perturbative QCD framework (^] . 

' An interesting question is whether the experiments at high-energy hadron colliders, such as the Tevatron or LHC, 
, can test saturation while for the moment this search is mainly considered for heavy-ion collisions. In the present paper, 
our aim is to look for saturation effects in the context of Mueller-Navelet lOj forward-jet production in hadron-induced 
, hard collisions. The key difference with electron-induced reactions is that the hard probe is no more a virtual photon 

7* but a virtual gluon g* , see Fig. la. 
\ A basic ingredient of the GBW saturation models is the QCD dipole formalism 0, ^| in which the hard cross- 
^ • sections read 

(N : 

' where rj=i^2 are the transverse sizes of the dipoles and A77 is the pseudo-rapidity range of the dipole-dipole cross- 
■ section ^^^(Ary, ri, r2). In our notations, (^^'-"(fijQi) are the dipole distributions in the target and projectile, and Qi 
' the virtualities of the hard probes that set the perturbative scale. 

Formula expresses a factorization property which has been shown to be equivalent Il3!| to fc^-factorization Ti'l 
in the BFKL framework. In this framework, the distributions (j)'-'^\ri,Qf) are related to the "impact factors" which 
Qh. describe the coupling of the target and projectile to the BFKL kernel. In the case of 7*— induced reactions, the dipole 
^ ' distribution functions 0^(r, Q^) are well-known from QED and the equivalence with photon impact factors checked. 

In the case of forward-jet production with transverse momentum > Q ^ 1 GeV, the corresponding distribution 
!>■ I (f>{r, Q^) can be derived |15llla| in the coUinear approximation i.e. in the Double Leading Log approximation (DLL). 

In the present paper, following an approach jl7l | for 7* — 7* cross-sections, we shall describe the predictions of 
saturation for g*— induced reactions such as Mueller-Navelet forward-jet production. For this sake, we will make 
' use of the dipole distribution <j){r,Q^) derived in [TsI ITgI . For simplicity, we assume the validity of the same GBW 
cross-sections as for 7*— 7*. The /cy-factorization property is assumed to be preserved in the presence of saturation (see 
a recent discussion in |l3|)- The question of going beyond this simple scheme e.g. using a more complete formulation 
of saturation is left for further work. 

The plan of the paper is as follows. In section II, recalling the results of [TEl ll^ . we show how the emission 
of a forward gluon jet can be recast in terms of a dipole distribution. We also present the GBW formulation of 
the dipole-dipole cross-sections. In section HI, we derive our results for the Mueller-Navelet jets cross-sections with 
saturation. In section IV, we discuss these results in the prospect of experiments at the Tevatron and LHC and 
propose characteristic observables for saturation. The final section V is devoted to conclusion and outlook. 
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FIG. 1; Forward jets and impact factors at hadron colliders. Fig la: Mueller-Navelet jets at hadron colliders. Fig lb: onium 
(+jet) impact factor in the partonic representation. Fig Ic: onium (+jet) impact factor in the dipole representation. A77 : 
rapidity gap between the two gluon jets, qr > Qi,2 '■ tranverse momenta of the gluon jets, k, k' : transverse momenta of 
the gluons interacting with the BFKL kernel. The gluon-dipole coupling {f°{k'^,r), see text) is sketched by the black point in 
Fig.lc. 



II. DIPOLE FORMULATION 
A. Forward jets and dipole distributions 

Let us first recall how one can obtain the dipole distribution 0(r, Q^) associated with a forward jet with transverse 
momentum qx > Q- The derivation is made |16j using the example of a final-state gluon being emitted from an onium 
{qq state) of size tq, see Figlb. QCD factorization will allow to extend the result to the case of an incident hadron. 

Assuming the condition 1 GeV^^ ^ tq ^ the onium is small enough to allow a perturbative QCD calculation 
but large enough with respect to the inverse transverse momentum of the forward jet. Using fc^-factorization in the 
BFKL framework the (unintegrated) gluon density /(fc^,ro) entering at each vertex of the BFKL cross-section (see 
Fig. lb) can be factorized in the following way: 

/(fc2,ro)^alogi / ^e{<f~Q^)f{\k + <^\^T^)^ T (2alogilog^l .h{Qr) f{k\r) (2) 

in the coUinear approximation rofj'T S> 1 for the onium. k is the transverse momentum of the gluon connected to the 
BFKL kernel (see Fig.l). /(/c^,ro) is the lowest order BFKL equation written in an unfolded form (see for instance 
[20j |) and the initial gluon density reads f^{k^,r) = 2a{l — Jo{kr))/k^ . Jq and Ji are the Bessel functions. 

Equation (O can be interpreted as the extension to forward jets of the equivalence between the momentum-space 
(partonic) and coordinate-space (dipole) representations. The middle term corresponds to the contribution displayed 
in Fig. lb. The last term is described in Fig.lc and matches with the Mueller picture of cascading dipoles in the 
1/Nc limit, in which the QCD wave function of an initial onium is expanded over multi colorless-dipole configurations. 
The factor in brackets { 2a log ^ log ^ } corresponds to the first order contribution of the Dokshitzer-Gribov-Lipatov- 
Altarelli-Parisi (DGLAP) gluon ladder 9], i.e. the probability of finding a dipole of size r inside the onium of size tq, 
at the Double Leading Log (DLL) approximation; thanks to QCD factorization properties, it is included in the gluon 
structure function of the incident particule. /° {k"^ , r) is nothing else than the factorized gluon density 0, |^ inside 
the dipole of size r which, in the dipole formulation is included in the dipole-dipole cross-section. 

Having factorized out both the contribution to the structure function and the one to the dipole-dipole cross-section, 
one is left with the function (/)(r, Q^) which describes the resulting size distribution of the interacting dipole. Hence, 



one is led to identify^ 



cj,{r,Q^) = f- MQr) . (3) 

The forward-jet emission is thus put in correspondence with a small colorless dipole of size r = 0{1/Q). The 
distribution of sizes around that value is given by 0(r, Q^) in (|^. 

Some comments are in order. Via its description in terms of dipoles, fcT-factorization leads to a description of the 
forward jet (coming from a colorful virtual gluon g*) in terms of colorless qq dipoles. Indeed, within the 1/Nc scheme 
of the dipole formalism, the color neutralization of the forward jet is described by a cascade of dipoles, as pictured in 
Fig.lc. This means that, in this representation, the color quantum number carried by the incoming virtual gluon g* 
is neutrahzed through the cascade of dipoles. 

The obtained dipole distribution (/)(r, Q^) is not everywhere positive, the Bessel function oscillating in sign for 
rQ > 4, and therefore cannot be interpreted as a probability distribution. We interpret this feature as a breakdown 
of the coUinear approximation. Hence, in our framework, we have to check the positivity of the cross-sections, as will 
be discussed later on. For the Mueller-Navelet BFKL cross-section, positivity is satisfied by construction. 

B. Dipole-dipole cross-sections with saturation 

Let us recall the formulation of the GBW saturation model for dipole collisions. Initially, the GBW approach is 
a model for the dipole-proton cross-section which includes the saturation damping of large-dipole configurations. For 
the description of 7* — 7* cross-sections at LEP, see Ref. it has been extended^ to dipole-dipole cross-sections. 

The same saturation scale is considered for dipole-dipole and dipole-proton cross-sections. 

The parametrisation of this dipole-dipole cross-section is 

...(A.,rr,r.)^.o{l-exp(-^^)}, (4) 

where i?o(A7]) = e^^^'^''^'^''^^^ /Qq is the rapidity-dependent saturation radius and the dipole-dipole effective radius 
r'^g{ri,r2) is defined in such a way to satisfy color transparency, namely add fx r|^i 2 when — > 0. As in 
three scenarios for roff(ri,r2) will a priori be considered: 

l.r!, = -fM, 2.r!, . nMrl,rl) 3..,^, = min(.?, r^^) {l-f In ^g^} . (5) 

All three parametrisations exhibit color transparency. Cases 1 and 2 reduce to the original GBW model when one of 
the dipoles is much larger than the other and the model 3 corresponds to the dipole-dipole cross-section mediated by 
a two-gluon exchange jl^l- For the saturation radius Ro{Ari) we adopt the same set of parameters"^ as those in [a.ll7|. 



III. HARD CROSS-SECTIONS 



Let us derive the general formulae we get for the cross-sections Defining u — r2/ri and r'^j.j — r\fi{u), the 
three scenarios considered in (0) can be rewritten 

f I \ ""^ f I \ / if u < 1 f I / - logu) if u < 1 , . 

^'^""^^iTu^ /2(") = |i if^>l /3(«) = |i^l^g^ . (6) 

Then inserting 121 in formula leads to 

— = I-Q1Q2 / du / rdr Ji{rQi)Ji{ruQ2) e 
c^o Jo Jo 

Jo V 



(7) 



^ This formula can also be obtained Il5l for (l>(r, z,Q^), taking into account the energy fraction 2 shared between the quark and the 
antiquark of the dipole. However, since the dipole-dipole cross-sections we will consider are 2— independent, we only have to consider 
the distributions integrated over z. 

^ For our purpose, wo shall only use a high-energy approximation of the expressions quoted by the authors |17|. 

3 A = .288, A»yo = 8.1 for Qo = l GeV. 



after integration over r. Ii is the modified Bessel function of the first kind. Formula |7J gives the theoretical cross- 
sections within the GBW model for hard hadronic probes. 

Let us discuss our results. The dipole distribution 0(r, Q^) is not everywhere positive but this is not a priori an 
obstacle as long as the corresponding total cross-sections ^ stay positive. This is for instance realized by the BFKL 
cross-section |l5| . It is compulsory to verify whether or not this positivity is altered by saturation. By numerical 
inspection of formulae lO, we checked that the positivity constraint is verified. Qualitatively, this is due to the fact 
that the negative values of the Bessel functions in (O are present for large dipole sizes whose contributions are strongly 
reduced by saturation. 

Another constraint is to check that the cross-sections add ^ '^cff /4^o(^'7) corresponding to the limit of small 
dipole sizes in lead to cross-sections behaving like 1/ {i?Q(A77) max (Q^, Qj)} -i expected from transparency. 
Computing the gluon-gluon cross-section in this limit gives 

o-o "^QlQl CTo ,2 n2^ c^o . / 1 1 \ 



which shows that the models 1 and 3 of (|SJ| verify the constraint. The model 2 does not, as confirmed by an explicit 
integration of lO which gives in this case 

at large Rq- Hence, within our approximations, the model 2 cannot be considered. 

It is possible to derive a general formula for saturation in the dipole framework which could be valid for any hard 
probe expressed in terms of the dipole basis, be it a forward jet, an onium, a virtual photon, etc... In particular, it 
will be useful to extend the saturation discussion to forward jets at HERA, with a 7* — g* cross-section. We consider 
the Mellin transforms of the dipole distributions 0(t) — J (Pr [r'^Q'^ Y (t){r,Q'^) . For instance, one has 

' r(l-r) ^ ^ ' r(4-2r)r(2 + 2T) ^ ' 

After some straightforward algebra, we obtain the following inverse Mellin transform expressions: 

- = / :^0(i)(r)(2Qii?o)-'" / ;^^(2)(a)(2g2i?o)^'" g,(a,r) ; < Re{a), Re{T), Re{a + r) < 1 , (11) 

ao J 2,111 J 2i'K 

where ^^^^ and 0^^-' are the Mellin-transformed dipole distributions in the target and projectile and, for the different 
models (01, one has 

r(l-T-a)r(a)r(r) r(l-r-a) 

^^^"'") = r(l + r + a) ^^^"'") = Vr (12) 

g3(fT,T) = -2— -r(-r-a){e2V-i-^-'^r(T + a + l,2a) + [r ^ a]} . 

This formulation allows us to extend easily our computations to various cases. After easy transformations, one gets 
for the different GBW models 

— = dx Ji(x) A(i)(xQii?o)^(2HxQ2i?o) ; A{x) = / -^x-^^~^{t)T{t) 

CTo Jo J 2«7r 

^ = r2xdx e-^" B^^\2xQiRo)B^^\2xQ2Ro) ; B{x) = / ^^-^-^ (13) 

ao Jo J 2«7r r 

2xdx e-^' (^C^^'>{2xQiRo)D(^\2Q2Ro,x)+D^^\2QiRo,x)C^^\2xQ2Ro)) 



CT3 _ 
Co o"o 



where 



CW./|L.-^M; D(Q..)^I^JQ.r-^j0^y (14) 

With these general formulae, we find back our previous results ()7I9|I and obtain* those for 7*. 



* Note that the input functions A, B, C, D I13I14I correspond to specific Meijer functions. 




FIG. 2: g* — g* cross-sections (model 3). Qi — Q2 = Q '■ symmetric virtuality case. Arj = (from left to right) 4,6,8,10 : 
rapidity intervals. Full lines: saturation cross-sections Q. Dashed lines: without saturation 



IV. PHENOMENOLOGICAL APPLICATIONS 



Let us investigate the phenomenological outcome, for hadron colliders, of our extension of the GBW models to 
hadronic {i.e. g*) probes. The theoretical cross-sections are obtained from formulae I©-© and (|11|I - H14I) . in terms 
of the physical variables Qi, Q2 and A77, once the saturation scale parameters Qq, A and Ar^o are taken identical to 
their reference values (see footnote 3). 

In Fig. 2, we display the cross-section ratio 0-3/ ctq as a function of \og{Q'^/QQ), where Q = Qi — Q2 for values of 
the rapidity interval Ary = 4,6,8, 10. We also compare with the corresponding ratio for the non-saturated case 
As expected, the curves show the well-known trend of the GBW model, namely a suppression of the non-saturated 
cross-sections, with a convergence towards the full saturation limit a ^ a^.lu order to appreciate more quantitatively 
the influence of saturation, it is most convenient to consider the quantities TZi/j defined as 

_ a(Qi,Q2,A77,) 
^/■' -^(Qi,Q2,Ary^.) ' ^^'^ 

i.e. the cross-section ratios for two different values of the rapidity interval. These ratios display in a clear way 
the saturation effects. They also correspond to possible experimental observables since they can be obtained from 
measurements at fixed values of the virtual gluon light-cone momentum and thus are independent of the gluon 
structure functions of the incident hadrons. Indeed, such observables have been used for a study of Mueller-Navelet 
jets for testing BFKL predictions at the Tevatron [Sill 111. 

In Fig. 3 we plot the values of 7?.4.6/2.4 (resp. TZgu) as a function of Qi = Q2 = Q- These ratios correspond to values 
for Mueller-Navelet jets studied at the Tevatron [22. (resp. realistic for the LHC '2^). The results are displayed 
both for models 1 and 3, see (0). 

As expected from the larger rapidity range, the decrease of TZ between the transparency regime and the saturated 
one is larger for the LHC than for the Tevatron. The striking feature of Fig. 3 is that the effect of saturation appears 
as a sharp transition for some critical range in Q (higher for the LHC). 

Let us compare the resulting ratios for hadronic probes ((7*-initiated) to those for the virtual photon (7*-initiated) for 
the same values of the rapidity ranges, see Fig. 3. Interestingly enough, the photon transition curve is much smoother, 
a phenomenon which can be explained by the different structure of the dipole distribution function. Indeed, as 
discussed previously (f^, the dipole distribution <j){r,Q^) has a tail extending towards large dipole sizes, which are 
more damped by the saturation corrections. Hence ip^r, Q^) is more abruptly cut by saturation than the photon dipole 
distribution <j)'^{r, Q^). 

In Fig. 4, we display the variation of TZs/4, when one looks for asymmetric situations, i.e. Qi > Q2- As seen from 
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FIG. 3: Cross-section ratios TZi/j- The resulting ratios for models 1 and 3 are plotted for rapidity intervals i = 8,j =4 and 
i = 4.6, j = 2.4. The comparison is made with 7*— 7* ratios for model 1 and equivalent kinematics. The non-saturated case 
would be a constant corresponding to the high limit of the plots. 



the figure, the transition may become even sharper in this case, with the formation of a bump at a rather high value 
of Qi, which could provide an interesting signal for the saturation scale. The origin of this bump lies in the different 
rate of increase of the cross-sections towards saturation when the virtualities are different. This possible signal is 
present at rather high scale, which could be useful for experimental considerations, as we shall develop now. 

i) Mueller-Navelet jets 

Two jets separated by a large rapidity interval, or Mueller-Navelet jets the more natural process for our 

formulae lO to be applied. Indeed, a measurement of those dijet cross-sections has been performed at Tevatron, with 
jets of transverse momentum with a lower Et cut, related to the virtuality Q. To actually mesure the ratio TZ in 
order to get rid of uncertainties on the structure functions, the two available incident energies (630 GeV and 1800 
GeV) were used. The result was a strong increase of TZ with the rapidity interval, which was pointed out as a possible 
hint of BFKL evolution. Saturation studies are favored by large rapidity intervals (as demonstrated in Fig. 3). The 
relevant range of virtuality Q for expecting a clear saturation signal is albeit rather low (see Fig. 3, 4). 

If the strong experimental signal reported in [22:] appears to be confirmed, the saturation prediction displayed in 
Fig. 3 could well be relevant (with a redefinition of the parameters). However, the BFKL evolution itself at Tevatron 
ener gies appears to be quite sizeably modified by finite energy, running coupling corrections p^l ] and by experimental 
cuts (2^1 . Anyway, a simulation of Mueller-Navelet jets that would incorporate the relation between the Et cut of 
the jet and the virtuality Q is needed to discuss the feasability of saturation tests in this case. 

ii) Heavy vector mesons 

As an alternative to hard forward jets, one could consider '27] the detection of two heavy vector mesons with 
moderate transverse momentum and separated by large rapidity intervals. Indeed, using J/'^'s or T's may provide a 
hadronic probe of precise mass and transverse momentum. It potentially realizes a colorless qq probe and thus could 
give an information on the differential distribution of dipoles (j){r,Q'^), for instance on the dipole-size distribution. 
Moreover, the leptonic decays may facilitate the event selection. 

iii) Charmed and beauty hadrons 

A forward-jet detection corresponds to one of the qq partners of a dipole. The detection of a heavy flavored meson 
would give a similar interesting signal. One thus would look for the detection of heavy flavored mesons separated by a 
large rapidity interval. In particular, the detection of a D* on one side and a i?-meson on the other side would realize 
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FIG. 4: 7^8/4 1 asymmetric case. The curves are drawn for model 3. Note the bump at Qi/Qo ~ 3. 



interesting asymmetric configurations''^ such as seen on Fig. 4. One could also play with their transverse momentum 
cuts to vary the g* virtualities. 

These possibilities of realizing hadronic probes of saturation certainly deserve more studies in the near future. 
Simulations of these processes at Tevatron and LHC energies will give a quantitative estimate of the potential of 
hadronic colliders to reveal new features of saturation. 



V. CONCLUSION AND OUTLOOK 

Let us briefly summarize the main results of our study. 

We started from an extension of the Golec-Biernat and Wiisthoff saturation model to hadronic collisions. For this 
sake, we used a QCD dipole formulation [isl Il6l | of hard hadronic probes based on a fc^^-factorization assumption; 
such probes as forward (Mueller-Navelet) jets, heavy mesons, heavy flavoured mesons are initiated by off-mass-shell 
gluons. Our results are: 

i) A derivation of saturation predictions for total hard cross-sections at hadron colliders, e.g. the Tevatron and 
LHC. 

ii) Observables which possess a potentially clear signal for saturation at high rapidity intervals and gluon virtualities 
around the expected saturation scale in dipole-dipole interactions, see Figs. (2-4). 

iii) The suggestion of using, besides the well-known Mueller-Navelet jets, heavy vector and/or heavy flavored mesons 
to measure hard cross-sections and their transition towards saturation. 

There are quite a few open issues for the present formalism: 

On a phenomenological ground, it indicates a way for simulations in the framework of hadron colliders using the 
QCD dipole formalism which had been so useful in the HERA context. These phenomenological studies will tell us 
whether and how saturation could be present and checked at the Tevatron and/or the LHC. Note also that forward 
jets at HERA, which are initiated by g* — 7* configurations with large rapidity separation, could be interesting to 
investigate. Beyond the scope of the present paper, determinant experimental issues, like fighting against pile-up 
events, background studies, possibility of a direct access to the hard cross-section ratios Tli/j, etc... deserve to be 
explored. 



^ Note that asymmetric configurations arc preferable to avoid eventual non-BFKL logarithmic corrections |2f^. 



On the theoretical ground, a study going beyond the single qq basis for the hard probe dipole distribution is deserved. 
In particular, adding the gqq or few-dipole configurations is important to discuss the /cT-fa-ctorization assumption. A 
more complete study of saturation effects in the emission of an energetic gluon ,1^ deserves further investigation. It 
may also allow us to extend our formalism to diffraction processes, since such configurations appeared important for 
the GBW model 's^ at HERA. It is also possible to extend the theoretical analysis beyond the GBW formulation and 
to directly introduce solutions of the non-linear QCD evolution equations, which would have a BFKL (and not 
merely transparency) limit at low density. 
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